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Triple gauge boson couplings are measured from W-pair events recorded by the OPAL detector 
; at LEP at centre-of-mass energies of 183 - 209 GeV with a total integrated luminosity of 
680 pb -1 . Only CP-conserving couplings are considered and SU(2)xU(l) relations are used, 
resulting in four independent couplings, « 7 , g\, A 7 and g\. Determining each coupling in 
a separate fit, assuming the other couplings to take their Standard Model values, we obtain 
K 7 =0.88i[J;^, £i=0.987lH!t, A 7 =-0.060±[X and ^=— 0.04±g;l|, where the errors include both 
statistical and systematic uncertainties. Fits are also performed allowing some of the couplings 
c5 ' to vary simultaneously. All results are consistent with the Standard Model predictions. 
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1 Introduction 



The triple gauge boson vertices, WW7 and WWZ, are expected in the Standard Model due to its 
non-Abelian nature. The Standard Model predicts exact values of the associated parameters, 
the triple gauge boson couplings (TGCs). Any significant deviation of the measured values 
from this prediction would be a signature of new physics beyond the Standard Model. Several 
measurements of these couplings have already been performed at LEP using W-pair [1-5] , single 
W [6] and single photon [7] production. Limits on TGCs also exist from studies of di-boson 
production at the Tevatron [8]. This paper reports results on TGCs using W-pair events from 
the full OPAL LEP2 data sample. 

Assuming a general Lorentz-invariant Lagrangian which satisfies electromagnetic gauge invari- 
ance, charge conjugation (C) and parity (P) invariance, five independent TGC parameters are 
expected to describe the WW7 and WWZ vertices. These can be taken A z and 

A 7 [9-12]. In the Standard Model gf=K, z =K,^=l and A z =A 7 =0. When the restriction on C and 
P invariance is dropped, but their product, CP, is required to be invariant, a sixth coupling, 
g\ [9-12], which violates both C and P, is added. This coupling vanishes in the Standard Model. 
Couplings which violate CP [13] are not considered in this analysis. 

Precision measurements on the Z° resonance and lower energy data motivate the following 
SU(2)xU(l) relations between the five C and P conserving couplings [9,11], 

Ak z = -A/t 7 tan 2 9 W + Agl, 
A z = A 7 . 

Here A indicates a deviation of the respective quantity from its Standard Model value and 9 W 
is the weak mixing angle at tree level, defined by cos9 w = M-w/Mz- These two relations leave 
three independent C and P conserving couplings, A/t 7 , Agl and A(=A 7 =A Z ), plus the C and P 
violating coupling g\. These couplings are not significantly restricted [14,15] by measurements 
from data collected at the Z° pole at LEP and SLC. 

While the single W and single photon processes can also be used to measure TGCs, W-pair 
production is by far the most sensitive process for this measurement. Its sensitivity comes 
from the triple gauge boson vertex connecting the intermediate s-channel Z°/7 to the outgoing 
W-bosons. The TGC contribution depends on the helicity states of the outgoing W bosons, 
determining the angular distributions of the W bosons and of their decay products. The total 
W-pair cross-section is also affected by the TGCs, but for centre-of-mass energies well above 
the threshold for W-pair production and for small values of anomalous couplings its sensitivity, 
compared with the angular distributions, is much lower. Very often, however, the angular 
distribution analysis yields two ambiguous solutions for the TGCs, one being quite far away 
from the Standard Model values [16], and the information from the total cross-section can help 
in resolving this ambiguity. 

The production and decay of W-pair events can be characterized by five angular variables. 
These are the W~ production polar angle 1 , 9w, the polar and azimuthal angles, 9* and 0*, 

1 The OPAL right-handed coordinate system is defined such that the origin is at the centre of the detector, 
the z-axis is parallel to, and has positive sense along, the e~ beam direction, 9 is the polar angle with respect 
to z and <p is tne azimuthal angle around z with respect to the x-axis, which points to the centre of the LEP 
ring. 
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of the decay fermion from the W~ in the W~ rest frame 2 , and the corresponding polar and 
azimuthal angles of the decay anti- fermion from the W + . The experimental accessibility of 
these angles depends on the final state produced when the W bosons decay. 

All W-pair final states are used in this study, namely the leptonic, tv[ivi>i, the semileptonic, 
qqive, and the hadronic, qqqq final states, with branching fractions of 10.5%, 43.9% and 45.6% 
respectively in the Standard Model. 

The next section discusses the OPAL data and Monte Carlo samples. The following three 
sections then present the first part of this study using the W-pair angular distributions. For 
this purpose, the kinematic observables of each event are reconstructed as described in Section 3 
and then used in Section 4 to construct optimal observables and analyse them in terms of the 
TGCs. The study of systematic errors for this part of the analysis is summarised in Section 5. 
Section 6 describes the second part of this analysis, where the total W-pair event rate is used 
to extract additional information on the TGCs. The TGC results are presented in Section 7 
and summarised in the last section. 



2 Data and Monte Carlo samples 

The OPAL detector is described in detail elsewhere [17]. The data were collected during 1997- 
2000 around eight different centre-of-mass energies. The integrated luminosity at each energy is 
evaluated using small angle Bhabha scattering events observed in the silicon tungsten forward 
calorimeter. The luminosity-weighted average values of the centre-of-mass energies and the 
corresponding luminosities are listed in Table 1. 

The main Monte Carlo generator used to simulate the W-pair signal is KandY, which is 
a combination of KORALW1.51 [18] and YFSWW3 [19] running concurrently. These pro- 
grams generate all four-fermion final states using the full set of electroweak diagrams includ- 
ing the W + W~ production diagrams (class 3 CC03) and other four-fermion graphs, such as 
e + e~ — ► Wez7 e , e + e~ — > Z°e + e~ and e + e~ — > Z°Z°. All diagrams are corrected for initial and 
final state radiation, while a more complete correction of 0(a) radiative effects is applied to 
the CC03 diagrams. For each centre-of-mass energy, 500,000 events have been generated us- 
ing these Monte Carlo programs, and these events are used as our reference samples for the 
TGC analysis. The TGC dependence is obtained by reweighting these samples using the four- 
fermion matrix element calculation taken from the EXCALIBUR [20] Monte Carlo program. 
This program has also been used to generate samples with anomalous TGCs which are used to 
cross-check our TGC extraction method. Although EXCALIBUR includes less complete 0(a) 
radiative effects, it can still be used for small anomalous TGC values, where the effect of the 
missing 0{a) corrections on the anomalous TGC contribution can be neglected. 

Background sources to the W-pair signal are four-fermion final states, such as qq££, which are 

2 The axes of the right-handed coordinate system in the W rest-frame are defined such that the z-axis is along 
the parent W flight direction in the overall centre-of-mass system, and the y-axis is in the direction e"xf 
where e~ is the electron beam direction and W is the parent W flight direction. 

3 In this paper, the three lowest order W-pair production diagrams, i.e. i-channel v c exchange and s-channel 
Z°/7 exchange, are referred to as "CC03", following the notation of [9]. 
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produced only by diagrams that do not involve the triple gauge boson vertex, two-fermion final 
states and two-photon processes. The four-fermion background is generated using KORALW1.42 
[21] and grc4f [22]. As an alternative for systematic studies we replace KORALW1.42 by 
EXCALIBUR. Two-fermion final states are generated using KK2F [23]. For systematic 
studies, KORALZ [24], for e + e~-*fi + fj,-, e+e'-^r+r- and e+e'^uu, and BHWIDE [25], 
for e + e~— >e + e~, are also used. Background to the semileptonic final state from single-tag 
two-photon processes, where the outgoing electron or positron is detected, is evaluated using 
HERWIG [26], and F2GEN [27] is used as an alternative for systematic studies. Background 
from leptonic final states from untagged two-photon processes, where both outgoing electron 
and positron escape undetected into the beam pipe, is calculated using the Vermaseren [28] 
and BDK [29] generators. 

To estimate the fragmentation and hadronisation systematics, Monte Carlo samples of W + W~ 
events were produced by the KORALW1.42 generator, with the fragmentation and hadronisa- 
tion stages generated separately by either PYTHIA5.7 [30], HERWIG6.2 [26] or ARIADNE4.08 [31]. 
These hadronic simulation generators have been tuned to Z° hadronic decays [32]. Similarly, for 
the Z°/7 — > qq background, the same hadronic simulation programs have been used to fragment 
and hadronise events produced by the KK2F generator. 

All Monte Carlo samples mentioned above were processed by the full OPAL simulation pro- 
gram [33] and then reconstructed in the same way as the data. 



3 W-pair event selection and reconstruction 

W-pair events are selected with the same procedure used in the W-pair cross-section measure- 
ment as described in references [3,34]. There are three independent selections corresponding 
to the three final states. Each candidate is then reconstructed in order to obtain the maximum 
possible information on the W production and decay angles, which are needed to extract the 
couplings. Events which cannot be well reconstructed are rejected from the sample, as will be 
described below. The numbers of candidates left at the different centre-of-mass energies for the 
three final states are listed, along with the expected values, in Table 1. 

Three kinematic fits with different sets of requirements are used in the event reconstruction: 

A. Require conservation of energy and momentum, neglecting Initial State Radiation (ISR). 

B. Additionally constrain the reconstructed masses of the two W-bosons to be equal. 

C. Additionally constrain each reconstructed W mass to the average measured value from 
the Tevatron 4 , M w =80.456 GeV/c 2 [35]. 

For qqqq events, where all four final state fermions are measurable, fits A, B and C have 4, 5 
and 6 constraints respectively. For qqeI7 e and qq/iX 7 ^ events the number of constraints is reduced 

4 The LEP results for the W mass are not used for the TGC measurement, since they have been obtained 
under the assumption that W pairs are produced according to the Standard Model, whereas W production at 
the Tevatron does not involve the triple gauge boson vertex. 
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JCdt 


qqive events 


qqqq events 


Ivil ' Vy 


events 


(GeV) 


(pb- 1 ) 


Observed 


Expected 


Observed 


Expected 


Observed 


Expected 


182.7 


57.4 


328 


331.1 


408 


418.9 


32 


37.4 


188.6 


183.0 


1090 


1123.5 


1437 


1388.3 


130 


124.3 


191.6 


29.3 


168 


182.8 


223 


222.8 


19 


19.9 


195.5 


76.4 


513 


478.9 


637 


594.3 


55 


52.5 


199.5 


76.6 


451 


479.7 


557 


585.0 


52 


52.1 


201.6 


37.7 


230 


236.3 


296 


286.1 


32 


25.3 


204.9 


81.6 


475 


512.4 


578 


606.5 


47 


53.2 


206.6 


136.5 


899 


854.0 


1051 


1012.3 


92 


90.0 


all 


678.5 


4154 


4198.7 


5187 


5114.2 


459 


454.7 



Table 1: Observed and expected numbers of data candidates selected for the angular distribu- 
tion analysis after all cuts in the different final states and for different centre-of-mass energies. 



by 3 due to the invisible neutrino. For qqrz7 r events there is at least one additional unobserved 
neutrino from the r decay, resulting in a loss of one constraint. The momentum sum of the 
track(s) assigned to the r can still be used as an approximation of the r flight direction, relying 
on its high boost, but the r energy is left unknown. Finally for tU^l events, where none of 
the leptons is a r, there are two invisible neutrinos. Hence, six constraints are lost and only 
requirement C is applied. 

In the following we discuss the reconstruction of each final state separately. 
3.1 Reconstruction of qqiUi final states 

Starting from the sample used for the W-pair cross-section measurement, we select candidate 
qqe!7 c and qq/il^ events with a reconstructed lepton track. The track charge is needed to 
reconstruct the W~ polar angle. For the qqrz/ r events, either one track or a narrow jet consisting 
of three tracks is assigned as the r decay product. 

The OPAL tracking detectors are used to reconstruct the muon momentum and the electron 
direction, whereas the electromagnetic calorimeters are used to give a more accurate measure- 
ment of the electron energy. As explained above, the direction of r candidates is given directly 
by the r decay products, whilst the energy of the r is estimated using a kinematic fit. 

The remaining tracks and calorimeter clusters in the event are grouped into two jets using the 
Durham k± algorithm [36]. The total energy and momentum of each of the jets are calculated 
using the charged track momenta and calorimeter energies and correcting for double counting 
employing the method described in [37]. 

The different kinematic fits described above are used to improve the resolution in the five 
kinematic variables used for the TGC analysis. This resolution is correlated with the kinematic 
fit probabilities which indicate the event reconstruction quality. Therefore we use the results of 
the kinematic fits to classify the events into nine quality classes which are used later in the TGC 
analysis. The exact definitions of these classes and the fraction of accepted events in each class 
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Class 


Final 

state 


%of 
events 


Requirements on fits A, B, C 


Kinem. variables 
taken from 


a) 


qqez7 c , qq/jjy^ 


30.4 


P c > 0.20 


fit C 


b) 


qqez7 e , qq>z7^ 


10.1 


0.05 < P c < 0.20 


fit C 


c) 


qqez7 e , qq/il/^ 


6.9 


0.01 < P c < 0.05 


fit C 


d) 


qqez7 e , qq>z7^ 


4.7 


0.001 < P c < 0.01 


fit A 


e) 


qqez7 e , qq>z7^ 


10.0 


P c < 0.001 


fit A 


f) 


qqeU e , qqfXu^ 


8.3 


C failed, P A > 0.001 


fit A 


g) 


qqez7 c ,qq/il7 M 


1.2 


C failed, A failed or P A < 0.001 


directly measured 


h) 


qqrV T 


21.7 


P B > 0.20 


fit B 


i) 


qqrV T 


6.8 


0.025 < P B < 0.20 


fit A 



Table 2: Definition of the nine quality classes of qqiise events. Pa, Pb and P c are the proba- 
bilities of kinematic fits A, B and C. 

are listed in Table 2. The event population in each class is well modelled by the Monte Carlo. 
To optimize the resolution in the kinematic variables used in the TGC analysis they are taken, 
according to the quality class, either from the directly measured values or from one of the three 
kinematic fits, as detailed in Table 2. The qqrV T events which fail kinematic fit B or pass it 
with fit probability below 0.025 are rejected. This cut suppresses those qqr77 T events which are 
correctly identified as belonging to this decay channel but where the r decay products are not 
identified correctly, leading to an incorrect estimate of the r flight direction or its charge. The 
fraction of such events in the qqrV T sample is thus reduced from 18% to 12%. 

Finally, the kinematic variables are used to calculate the optimal observables (see Section 4) 
which are used to extract the TGCs. Those events (typically 1% of the sample) with optimal 
observable values in the far tails of their distributions are interpreted as being badly measured 
and discarded. 

Performing the complete analysis on Monte Carlo samples generated with different anomalous 
couplings shows that our event selection does not introduce any bias on the TGC results. 

The purities of the final samples are about 96%, 98% and 91% for the qqez/ e , qq/-" 7 ^ and qqrz/ T 
final states respectively. The background consists mainly of TGC-independent four-fermion 
final states and Z°/7 — > qq events. Cross migration between the qqru T and each of the qqez7 e 
and qq/iX 7 ^ decay channels is at the level of 4-5% and is considered as signal. 

In the reconstruction of the qqfcvz events we obtain cos #w by summing the kinematically fitted 
four-momenta of the two jets. The decay angles of the leptonically decaying W are obtained 
from the charged lepton four-momentum, after boosting back to the parent W rest frame. For 
the hadronically decaying W it is not possible to distinguish between the jets of the quark and 
anti-quark. This ambiguity is taken into account in the analysis described below. 

In Figure 1 we show the distributions of the five angles corresponding to the qqlVz data sample 
and the expected distributions for A = ±0.5 and 0. These expected distributions are obtained 
from Standard Model Monte Carlo samples. To produce the A = ±0.5 distributions the events 
have been appropriately reweighted using the EXCALIBUR matrix element calculation. All 
Monte Carlo distributions are normalised to the number of events observed in the data. The 
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OPAL 



s 600 




Figure 1: Distributions of the kinematic variables cos^w, cos 9}, <fi}, cos#* t and 0* t , obtained 
from the qqip£ candidates. The solid points represent the data. The histograms show the 
Monte Carlo expectation of the Standard Model (solid line) and the cases of A = +0.5 and 
A = —0.5 (dashed and dotted lines respectively). The Monte Carlo histograms are normalised 
to the number of data events. The shaded histograms show the non-qqiVi background. In 
the case of W + ^lv decays the value of §\ is shifted by 180° in order to overlay W + and 
distributions in the same plot, which is valid so long as CP is conserved. 
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shaded histograms show the non-qq£z7 £ background and their normalisation relative to the 
signal is according to the Standard Model. For the cos#* et and 0* ct distributions the jet with 
< 0* ct < 180° is arbitrarily chosen as the quark jet from the decay of the W~, or the anti-quark 
jet from the decay of the W + . Sensitivity to A is apparent mainly for cos#w The contribution 
of cos 9}, (f)}, cos6>* t and 0* t to the overall sensitivity enters mainly through their correlations 
with cos#w 



3.2 Reconstruction of qqqq final states 

Using the Durham k± algorithm [36] each selected qqqq candidate is reconstructed as four jets 
whose energies are corrected for double counting of charged track momenta and calorimeter 
energies [37]. These four jets can be paired into W-bosons in three possible ways. To improve 
the resolution on the jet four-momenta we perform the five-constraint kinematic fit B for each 
possible jet pairing. In the following we consider only good jet pairings, namely those pairings 
with a successful fit B yielding a W-mass between 70 and 90 GeV. We accept only events with 
at least one good jet pairing. 

After this cut the overall efficiency is around 80%, varying slightly with the centre-of-mass 
energy, and the contamination from other W + W~ decays is below 0.3%. The major background 
contribution is from Z°/7 — » qq where additional quarks or gluons are radiated off the primary 
quarks. This contribution is between 12% to 17%, decreasing with centre-of-mass energy. The 
contribution of TGC-independent four-fermion final states is 4.5%, except at 183 GeV where it 
drops to 3%. 

We use a jet charge method as an estimator for the W charge. For the two jets coming from 
one W candidate we calculate the momentum weighted charge of all tracks contributing to this 
W candidate as, Q w = %|P||,i| ' 5 /^lf=i |P||,i| ' 5 , where iV w and N tot are the numbers of 

th 

tracks in the W candidate and in the total event respectively, qi is the charge of the i track 
and piu is the momentum component parallel to the jet axis. We take the W candidate with 
the more positive charge as the W + . 

Around 38% of the events have more than one good jet pairing, but the fraction of events where 
all three possible jet pairings are good is below 5%. To select the correct jet pairing, the good 
pairings in the event are processed through a likelihood algorithm [38] which takes into account 
correlations between the input variables. We use as input to our likelihood the W candidate 
mass as obtained from the kinematic fit B, the absolute value of the difference between the W 
candidate masses as obtained from the kinematic fit A and the absolute value of the jet charge 
difference between the two W candidates. The algorithm has been tuned on Monte Carlo events 
at each centre-of-mass energy in order to give the optimal performance. The di-jet pairing with 
the highest likelihood output is taken as the correct pairing. The distribution of the highest jet 
pairing likelihood output is shown in Figure 2(a). The distribution of the W charge separation 
for all highest likelihood pairings is shown in Figure 2(b). 

According to the W + W" KORALW Monte Carlo, the probability of selecting the correct di-jet 
combination is about 80%, and the probability of correct assignment of the W charge, once 
the correct pairing has been chosen, is about 77%. Both probabilities vary slightly with the 
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IQwi"Qw2' 

Figure 2: Distributions for qqqq data events (points) for the KandY four-fermion Monte Carlo 
(histogram) of (a) jet-pairing likelihood corresponding to the most likely combination; (b) 
charge separation between the two W candidates. 
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centre-of-mass energy but depend strongly on the jet pairing likelihood output. 

In order to reflect this dependence on the jet pairing likelihood output we use the output as 
a quality variable to classify the events into four quality classes which are later used in the 
TGC analysis. The definitions of these classes and the fraction of accepted events in each class 
are listed in Table 3. For classes with higher jet pairing likelihood output the probability for 
selecting the correct di-jet combination and the correct W charge assignment increases. The 
numbers are given in Table 3. The population of each class is well modelled by Monte Carlo. 



Class 


Jet pairing 


% of events 


Correct jet 


Correct charge 




likelihood output 




pairing probability 


assignment probability 


a) 


0.72 - 1.00 


27.2 


91.8% 


89.4% 


b) 


0.54 - 0.72 


24.9 


84.4% 


74.1% 


c) 


0.35 - 0.54 


22.4 


75.3% 


66.2% 


d) 


0.00 - 0.35 


25.5 


56.4% 


67.9% 



Table 3: Definition of the four quality classes of qqqq events. The correct charge assignment 
probability is calculated only for those events with correct jet pairing. 



The distributions of the variables cos#w, cos0* et and 0* t are shown in Figure 3, along with the 
expected distributions for the Standard Model (A=0) and for A=±0.5. In the cos0* et and 0* t 
distributions there are two entries per event, and according to our convention (see Section 3.1), 
the jet with < 0* et < 180° is chosen as the quark jet from the decay of the W~, or the 
anti-quark jet from the decay of the W + . The Monte Carlo histograms are normalised to 
the number of data events. The shaded histograms represent the separate contributions from 
wrong pairings, correct pairing but wrong charge and background. Their relative normalisation 
is according to the Standard Model. The data are again described well by the Standard Model 
prediction. 

3.3 Reconstruction of iuit'uir final states 

The reconstruction of ItUftv^ final states is possible, using kinematic fit C, only if there is no 
r-lepton in the final state. To reject events with r-leptons we use the lepton identification algo- 
rithm designed for the W branching ratio measurement [34]. The event rate analysis selection 
sometimes selects events with only one high momentum charged lepton reconstructed in the 
detector: such events were removed from the angular distribution analysis as the momenta of 
both charged leptons must be measured. As a result of these two cuts, the efficiency for tV(t 
events with £, £' — e, // drops from ~88% to ~73%, varying slightly with centre-of-mass energy. 
The contamination left in the sample from IVgl final states with a r-lepton is 14% and the 
background from other final states is 2%. 

Kinematic fit C has zero constraints and reduces to solving a quadratic equation. In the ideal 
case of no measurement errors and satisfying the conditions where fit C is valid, namely no 
ISR and narrow W width, one expects to obtain two real solutions corresponding to a two-fold 
ambiguity in the angle set cos 9w, <j>* an d 4>2- There is no ambiguity in the angles 6{ and Q\ which 
have a one-to-one correspondence with the lepton momenta. In realistic conditions, however, 
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Figure 3: Distribution of the kinematic variables cos^w, cos#* ct , and 0* ot for qqqq events. In 
the cos9* ct , and 0* t distributions there are two entries per event. The solid points represent 
the data. The histograms show the Monte Carlo expectation of the Standard Model (solid 
line) and the cases of X=+0.5 and A=-0.5 (dashed and dotted lines respectively). The Monte 
Carlo histograms are normalised to the number of data events. The open histograms show the 
contribution from correct pairing and correct charge assignment. The single hatched histograms 
represent the contribution from correct pairing, but wrong charge assignment. The cross- 
hatched histograms correspond to the contribution from wrong pairing, and the dark histograms 
show the contribution of background. 
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including ISR, finite W width and measurement errors, one may obtain no real solution, but a 
pair of complex conjugate solutions. The fraction of events with complex solutions falls with 
increasing centre-of-mass energy from 30% to 24%, in agreement with Monte Carlo prediction. 
In these events, the imaginary parts of the complex solutions are set to zero, yielding a single 
real angle set. The total weight of all events is the same, independent of whether one or two 
solutions are employed. In the latter case, the two solutions each are assigned a weight of one 
half. 

Figure 4 shows the distributions of the five reconstructed angles and the expected distributions 
for A=0, ±1. The Monte Carlo histograms are normalized to the number of data events. The 
shaded histograms represent the contribution of background and their normalisation relative to 
the signal is according to the Standard Model. The data distributions agree with the expectation 
of the Standard Model. 



4 Optimal observable analysis 

The angular distributions from all final states are consistent with the Standard Model expec- 
tation and no evidence is seen for any significant contributions from anomalous couplings. For 
a quantitative study we use the method of optimal observables in the same way as in [4]. This 
method relies on the linear dependence of the triple gauge vertex Lagrangian on the TGCs, 
corresponding to a second-order polynomial dependence of the differential cross-section, 

= S (0 \Q) + ]>> • SP(Q) +EW • ^(ft), a t = A« 7 , Arf, A and gl 



1,3 



where S<-°\ Sf > and S\f are functions of the five phase-space variables, ft=(cos#w! cos#*, <j)\, 
cos 6*2, 4>*£). For this kind of dependence it has been shown [39] that all information contained 
in the phase-space variables ft, which is relevant to the four couplings, is retained in the whole 
set of 14 observables Of ] = sf } (ft)/S(°)(ft), Of ] = (ft) / (ft) and Og } = Of = 

S\f (ft) / S^\VL) = 1,2,3,4). In this analysis we use the mean values of Of ] and { ff in 

addition to the mean values of o[ 1] [40]. 

In the one-parameter fit of ctj we only use the corresponding first order, and second 

(2) 

order, 0\ , observables. In this way, the other couplings are not used in the fit. Similarly, 
in two-parameter fits where two couplings are allowed to vary simultaneously, we use only the 
five observables corresponding to these couplings and their mutual correlations. When three 
couplings are allowed to vary, nine observables are used. 

The calculation of the mean optimal observables and their expected values is done separately 
for each centre-of-mass energy and for each of the qq£77^, qqqq and iV(l final states. The 
optimal observables are constructed for each event k with the set of phase-space variables 
ftfc using the analytic expression for the CC03 Born differential cross-section to calculate the 
values of S'^(ftfe), S^(ftfc), S'^(ftfe) and S^(ftfc). This calculation takes into account the 
reconstruction ambiguities for each particular final state. For the qqqq channel, the ambiguity in 
the W-charge determination is taken into account by weighting the , and Sff functions 
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Figure 4: Distributions of the kinematic variables cos 9^, cos 9\, cos 9* 2 and 4> 2 in the ^Vfl vgt 
analysis, where cosO^, (fil (cos 9%, §\) are the decay angles of the charged lepton (anti-lepton) 
from the W~ (W + J in the parent W rest frame. All events enter with a total weight of one, 
but in the case of events with two ambiguous solutions, each solution enters with a weight of 
0.5. The histograms, normalized to the number of data events, show the expectation of the 
Standard Model and the cases of X = ±1. The shaded histograms represent the contribution 
of the background. 
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for each charge hypothesis with the corresponding probability. This probability is determined 
from the Monte Carlo as a function of the estimated charge difference \Qwi ~ Qw 2 \ between 
the two W candidates in the event. 

For qqive and qqqq events the mean optimal observables are calculated using weights assigned 
to the event according to its quality class (see Tables 2 and 3). For qq£&e events, the weight for 
a particular class is inversely proportional to the variance of the optimal observable distribution 
of the events in that class. For qqqq events, the weights are calculated in a similar way, using 
the variances of the optimal observable resolution distributions. In this way, events where the 
quality of the reconstruction is poor obtain a lower weight. Using this weighted mean has been 
found, as expected, to enhance our sensitivity. 

The expected mean values of each observable as a function of the TGCs are calculated from 
four-fermion reference Monte Carlo samples generated according to the Standard Model. The 
EXCALIBUR matrix element calculation is used to reweight the Monte Carlo events to any 
TGC value required. The contribution of background events is also taken into account in the 
calculation of the expected mean values using corresponding Monte Carlo samples. 

A x 2 fit of the measured mean values O to the corresponding expectations E[0](a) is performed 
to extract the couplings from the data. The covariance matrix for the mean values is calculated 
from the Monte Carlo events and scaled to the number of data events. 

In some of our fits the x 2 function has a double minimum and therefore it is necessary to 
validate the errors obtained from the x 2 fits- This is performed with a large number of Monte 
Carlo "experiments". In each of these "experiments", we use for each centre-of-mass energy 
and each final state a subsample of Standard Model Monte Carlo events. The size of each 
subsample corresponds to the luminosity for that final state at the particular centre-of-mass 
energy. Background events are also included in these subsamples in the appropriate proportions. 
For each "experiment" , the various subsamples are analysed in the same way as the real data 
events. The distributions of the fit results are centred around the expected values, but there 
are some non-Gaussian tails. We therefore test the reliability of the error interval, as given by 
the region where x 2 is 110 more than 1 above its minimum value, by counting the fraction of 
"experiments" where the correct value falls within this region. The error estimate is considered 
to be reliable if the calculated fraction is consistent with 68%; otherwise the corresponding 
elements of the covariance matrix for the O values are scaled up by the necessary factor to 
obtain 68% of the subsamples within the error interval. The resulting scaled covariance matrix 
is also used to analyse the real data. This is done separately for each fit. The resulting scale 
factors vary between 1 and 1.076 for the one-parameter fits. However, in the three-parameter 
fit of Ak 7 , Agf and A, where double minima occur more frequently, a larger scale factor of 
1.28 is obtained. These scale factors have been found to be appropriate also for the 95% C.L. 
intervals. 



5 Systematic errors in the angular distribution analysis 

The following sources of systematic uncertainty are considered, as listed in Table 4. 
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There is a theoretical uncertainty in the 0(a) correction due to missing higher orders 
in the YFSWW3 Monte Carlo generator program. As a conservative estimate of this 
uncertainty, this O{o) correction is removed, degrading the events to KORALW1.42 level. 
The efffect of this change is taken as a systematic uncertainty. 

KORALW includes ISR up to 0(a 3 ). As a conservative estimate of missing higher orders, 
we replace it with ISR up to 0(a) only and take the difference as a systematic uncertainty. 

Another uncertainty in the signal Monte Carlo generator is related to the PYTHIA frag- 
mentation and hadronisation model which is used. To study the effect of this uncertainty 
we compare results of the analysis when applied to W + W~ Monte Carlo samples generated 
with KORALW1.42 using either PYTHIA, HERWIG or ARIADNE for the fragmentation 
and hadronisation phase. The effect of changing to HERWIG is the larger and is used 
to estimate the uncertainty. The main effect is on the qqqq final state, where HERWIG 
predicts a higher probability of correct jet pairing and correct W charge assignment than 
PYTHIA. This effect is partially due to the average charged multiplicity for light quarks 
being lower in our tuned version of HERWIG than in PYTHIA. Since the data are in 
agreement with PYTHIA [41], we weight the HERWIG Monte Carlo events so as to repro- 
duce the same charged multiplicity distribution as PYTHIA. The effect of this weighting 
is a reduction of the systematic error by about 15-30%. 

To assess the uncertainty in the TGC matrix element calculation which is used to weight 
our Monte Carlo events, we replace the calculation from the EXCALIBUR Monte Carlo 
generator program with the one from grc4f . 

The Monte Carlo simulation of the OPAL detector has been studied using back-to-back 
jets and leptons in Z° events collected during calibration runs which were taken regularly 
every year between the high energy runs. Some scaling and smearing had to be applied 
to the reconstructed jets and leptons in the Monte Carlo in order to achieve the best 
agreement with data. The uncertainties in these corrections correspond to variations in 
our results which are taken as systematic errors. These corrections are usually applied 
after the event selection. To estimate their effect on the selection of our sample, they 
were applied as a systematic check before the event selection, and the difference from 
the standard result was added in quadrature with the effects of the uncertainties on the 
corrections. 

Uncertainties in the background estimation are determined by varying both its shape 
and normalisation. To study the effect of the shape we use alternative samples for each 
background source, keeping the normalisation fixed. For the TGC-independent four- 
fermion background we use EXCALIBUR and grc4f instead of KORALW. For the 
Z°/7 -> qq background we use KK2F with HERWIG or ARIADNE rather than PYTHIA 
fragmentation. Finally, the HERWIG two-photon samples are replaced with F2GEN. For 
the background normalisation uncertainty the contribution from each source is scaled, one 
at a time, by the factors 1.2, 1.3 and 2.0 respectively. These scale factors are determined 
from studies with control data samples and from the differences between the alternative 
Monte Carlo generators, as explained in [34] . 

The consequences of the uncertainties in the centre-of-mass energy and the W mass are 
estimated as follows. For the centre-of-mass energy uncertainty of 40 MeV, Monte Carlo 
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samples at different centre-of-mass energies are used as pseudo-data. For the W mass we 
assume M w =80.456±0.059 GeV/c 2 [35] which is 0.126 GeV/c 2 higher than the one used 
in our Monte Carlo generators. We do not correct our TGC fit result for this shift in Mw, 
but it is accounted for in our systematic error, along with the Mw measurement error. 
The combined effect is assessed using Monte Carlo samples generated with different W 
masses. 

h) Bose-Einstein correlations (BEC) in qqqq events might affect the reconstruction of the W 
bosons and measured W charge distribution. We investigate this effect with a Monte Carlo 
program simulating BEC via re- adjustment of final state momenta using the LUBOEI 
model [42]. When both inter-W and intra-W correlations are taken into account, the bias 
is found to be larger than the case where only intra-W boson correlations are present. 

i) The jet reconstruction and the measured W charge distribution in qqqq events might also 
be affected by colour reconnection. This effect is investigated with several MC samples 
generated according to the Sjostrand-Khoze [43] models. We consider model I with colour 
reconnection probabilities of 34% and 97%. Models II and II' are also considered. Model 
I with colour reconnection probability of 97% produces the largest bias. 

The Monte Carlo samples are larger than the data by a factor of 150 - 900, depending on the 
centre-of-mass energy. Therefore, the uncertainty due to Monte Carlo statistics is neglected. 

The effect of each of these sources on the mean values of the optimal observables is estimated 
separately for each centre-of-mass energy and for each final state. The deviations obtained in 
the mean observables are used to construct a systematic error covariance matrix, where the 
systematic deviations from each source are assumed to be fully correlated between all centre- 
of-mass energies. For the combination of all final states, the covariance matrix is extended 
to incorporate the three final states. All uncertainties, except for those associated with the 
background, are assumed to be fully correlated between the relevant final states. This matrix is 
added to the statistical covariance matrix, and the fit is redone to obtain the TGC parameters 
with their total (statistical and systematic) errors. 



Source 


Ak 7 


A<7? 


A 


9l 


a) 


0(a) correction 


0.029 


0.011 


0.010 


0.031 


b) 


ISR 


0.006 


0.002 


0.003 


0.003 


c) 


Fragmentation 


0.038 


0.013 


0.018 


0.047 


d) 


TGC matrix element 


0.011 


0.005 


0.004 


0.005 


e) 


Detector simulation 


0.006 


0.003 


0.005 


0.013 


f) 


Background 


0.015 


0.004 


0.011 


0.012 


g) 


yfe, M W 


0.008 


0.004 


0.005 


0.017 


h) 


Bose-Einstein correlations 


0.001 


0.000 


0.004 


0.005 


i) 


Colour reconnection 


0.004 


0.003 


0.005 


0.000 


Total 


0.053 


0.019 


0.026 


0.061 



Table 4: Contributions to the systematic uncertainties in the determination of the TGC param- 
eters. The total systematic errors are obtained by combining all contributions at the optimal 
observable level and then calculating the effect on the TGC parameters. 
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For demonstration, we also calculate the separate contribution from each source of systematic 
error to the error on the TGC parameters, after combining the different centre-of-mass energies 
and final states. These are listed in Table 4. 



6 Event rate TGC analysis 

The event rate study, unlike the angular distribution analysis, requires a detailed investigation 
of the systematic uncertainties associated with the overall selection efficiency This investigation 
is part of the W-pair cross-section analysis and has to be done separately for each centre-of- 
mass energy and for each final state. So far, it has been completed for centre-of-mass energies 
of 183 GeV [3] and 189 GeV [34]. The TGC event rate analysis is therefore also restricted 
to these two centre-of-mass energies. In the case where the angular distribution analysis has 
two TGC solutions, one being far away from the Standard Model value, the results from this 
restricted event rate analysis would still be sufficient to reject that solution. The data from 
higher centre-of-mass energies would make a negligible additional contribution to our results. 
Compared with the angular distributions, the sensitivity of the event rates to TGCs around the 
Standard Model values is much smaller and its dependence on the overall luminosity is weaker. 

The event rate analysis has been described already in our previous publications [3,4]. The 
numbers of expected events are re-evaluated for this analysis using the KandY Monte Carlo 
samples to obtain the Standard Model expected number of events. 

The same selection, efficiency and background estimates as in the total cross-section analysis are 
used here. The numbers of selected events in each final state and for each of the two centre-of- 
mass energies are quoted in Table 5 along with the Monte Carlo expectations. The systematic 
errors on the expected values are due to Monte Carlo statistics, 0.5% uncertainty in the total 
cross-section, 0.2% uncertainty in the luminosity, uncertainties due to data/MC differences, 
tracking losses, detector occupancy and fragmentation. A detailed description of all these 
sources can be found in [34] . The expected numbers of signal events include contributions from 
all four-fermion final states which can be produced by diagrams involving the WW7 and WWZ 
triple gauge boson vertices. The other four-fermion final states are considered as background 
along with the two-fermion final states and two-photon processes. 





Luminosity 


Final 


Observed 


Expected events 


(GeV) 


(Pb^ 1 ) " 


state 


events 


Total 


Signal 


Background 


182.7 


57.21 ±0.25 




78 


77.5 ±3.1 


76.1 ±2.2 


1.4 ±2.2 






qq£V e 


361 


364 ±6 


337 ±6 


27 ±3 






qqqq 


438 


436 ± 10 


344 ±5 


92 ±9 


188.6 


183.05 ±0.40 




276 


287.2 ±5.0 


272.6 ±4.1 


14.6 ±2.9 






qq£v e 


1246 


1248 ± 16 


1171 ± 13 


76 ± 10 






qqqq 


1546 


1497 ±26 


1189 ± 14 


309 ±21 



Table 5: Observed and expected numbers of events in each W + W~ final state for the centre- 
of-mass energies of 182.7 GeV and 188.6 GeV. The separation between signal and background 
is explained in the text. All the quoted errors are systematic. 
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The total numbers of expected events are consistent within the statistical and systematic errors 
with the observed numbers. Therefore, there is no evidence for any significant contribution 
from anomalous couplings. A quantitative study of TGCs from these results is performed by 
comparing the numbers of observed events in each of the three event selection channels at each 
of the two centre-of-mass energies with the expected numbers of signal and background events, 
where the number of signal events is parametrised as a second-order polynomial in the TGCs. 
This parametrisation is based on the second-order polynomial dependence of the cross-section, 
as was used also in the optimal observable analysis (see Section 4). The polynomial coefficients 
in terms of event rate for each final state are calculated by weighting our Monte Carlo events to 
correspond to various TGC values using a four-fermion matrix element calculation procedure 
taken from the EXCALIBUR Monte Carlo generator program. 

The probability to observe the measured number of candidates, given the expected value, is cal- 
culated using a Gaussian distribution where the errors also include the systematic uncertainties. 
We assume full correlation between the systematic uncertainties for the three final states except 
for efficiency and background which are taken to be uncorrelated. All systematic uncertainties 
are assumed to be fully correlated between the two different centre-of-mass energies. 



7 Combined TGC results 

The results of the various fits are listed in Table 6. We present first the results of the angular 
distribution analysis without systematic errors obtained by performing one-parameter fits for 
each final state simultaneously for the eight centre-of-mass energies. For comparison, the aver- 
age fit errors of the Monte Carlo "experiments" are also quoted in Table 6, demonstrating the 
different sensitivities for the various final states and TGC parameters. In most cases, the fit 
errors are rather close to the expected ones, except for A/t 7 , where the fit errors for some final 
states are also very asymmetric due to the proximity of a second minimum in the y 2 function. 

Next we present the results of the angular distribution analysis with the systematic errors 
included as explained in Section 5. This modifies the relative contributions to the fit of the 
different optimal observables and different centre-of-mass energies, as well as the correlations 
between them. Therefore, not only the fit errors are affected, but also the fit results themselves. 
The corresponding logL curves 5 for the different channels are shown in Figure 5. The results 
from the three final states are combined, taking into account the common systematic errors as 
explained in Section 5, and the combined results are listed in Table 6. 

The results of the event rate analysis are presented by log L curves shown by the thin solid lines 
in Figure 5. They are combined with the results of the angular distribution analysis and the 
combined log L plots are shown by the thick solid lines. The correlation between the results of 
the two analyses, which is due to uncertainties in the background normalisation, centre-of-mass 
energy and W mass, is below 7%. Therefore this correlation is neglected in the combination. 
The numerical results of the combination are listed in Table 6 along with the corresponding 
X 2 /d.o.f. As expected, the event rate information has relatively little impact on the final results. 

5 Following the convention in LEP, we prefer to display our results by the logarithm of the likelihood function, 
logi, related to the % 2 function by, -logL= x 2 /2. The plotted function is -AlogL, which is obtained by 
subtracting from -logi its minimum value. 
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Including systematics 
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Including event rate 


n i Q+0.09 
— U.1Z_ 08 


U.UlO_ 33 


-o.oeo+Hll 


-0.0418:11 


Fit xVd.o.f. 


43.4/53 


39.4/53 


36.2/53 


42.0/53 


Standard Model x 2 /d.o.f. 


45.1/54 


39.5/54 


39.2/54 


42.1/54 


95% C.L. limits 


[-0.27, 0.07] 


[-0.077, 0.054] 


[-0.13, 0.01] 


[-0.28, 0.21] 



Table 6: Measured values of the TGC parameters from the angular distributions of each WW 
final state, using centre-of-mass energies of 183 - 209 GeV and the results after combining all 
final states. For the individual final states we also list the results before including the systematic 
errors and the expected statistical errors. We also list the results after combining with the W- 
pair event rate information from centre-of-mass energies of 183 and 189 GeV. The quality of 
the fits is demonstrated by their respective x 2 values and numbers of degrees of freedom. The 
compatibility with the Standard Model is shown by the \ 2 values obtained when the TGCs 
assume their Standard Model values. 
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Figure 5: Negative log-likelihood curves obtained from the angular distribution information 
of W-pair events taken at the centre-of-mass energies of 183 - 209 GeV, separately for the 
three different final states: qq£V(, (dotted lines), qqqq (dashed lines) and tu e l vi> (dash-dotted 
lines). The thin solid lines describe the contribution of the event rate information from all final 
states corresponding to centre-of-mass energies of 183 and 189 GeV. The systematic errors are 
included. The thick solid lines are obtained by combining all sources of information. The curves 
for each TGC parameter are obtained by setting the other three parameters to their Standard 
Model values. 
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The 95% confidence level intervals are also listed in Table 6. To show the compatibility with 
the Standard Model we also list in the table the x 2 /d.o.f. for Standard Model TGC values. 



To study correlations between the C- and P-conserving TGC parameters, Ak 7 , Ag\ and A, we 
perform fits where two parameters are allowed to vary while the third parameter and g\ assume 
their Standard Model values and are not used in the fit, as explained in Section 4. Similarly, we 
also perform a three-parameter fit where all three parameters are allowed to vary, leaving out 
only g\. Figure 6 shows the 95% C.L. contour plots obtained from the two-parameter fits and 
the corresponding two-dimensional projections of the three-parameter fit. The contour plots 
show a double minimum structure. Adding information from the single W channel, where the 
-logL function tends to have a steep rise with increasing Ak 7 [6], is expected to resolve this 
ambiguity. 

A x 2 test can be used to study the compatibility of our data with the Standard Model, in the 
same way as was done for the one parameter fits. Inserting the Standard Model values for the 
TGC parameters we obtain x 2 — 211.6 for 222 degrees of freedom. 



8 Summary 

Using a sample of W + W~ candidates collected at LEP at centre-of-mass energies of 183 - 
209 GeV, we measure the CP-conserving triple gauge boson couplings re 7 , g\, A and g\. For 
this measurement we use both the angular information of W-pair events and their total event 
rate. The results obtained are: 



where each parameter is determined from a single-parameter fit, while the other three param- 
eters assume their Standard Model values. These results supersede those from our previous 
publications [1-4]. They are all consistent with the Standard Model predictions of 1, 1, and 
respectively. This measurement constitutes strong evidence for the gauge structure of the 
Standard Model. 
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Figure 6: The two-dimensional correlation contours corresponding to — AlogL = 2.995 for 
different pairs of TGC parameters, as obtained from two-parameter (solid curves) and three- 
parameter (dashed curves) fits. The value 2.995 corresponds to a 95% C.L. region in two 
dimensions. In the case of the three-parameter fit, the curve is a projection onto the two- 
dimensional plane of the envelope of the three-dimensional — AlogL = 2.995 surface. The stars 
indicate the Standard Model expectations. 
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